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PDE-ODE SYSTEM 
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Abstract. We study the asymptotic behavior of solutions of one coupled PDE-ODE system 
arising in mathematical biology as a model for the development of a forest ecosystem. 
, In the case where the ODE-component of the system is monotone, we establish the existence 

of a smooth global attractor of finite Hausdorff and fractal dimension. 

The case of the non-monotone ODE-component is much more complicated. In particular, the 
set of equilibria becomes non-compact here and contains a huge number of essentially discontin- 
uous solutions. Nevertheless, we prove the stabilization of any trajectory to a single equilibrium 
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, if the coupling constant is small enough. 
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1. Introduction 
We study the following coupled ODE-PDE system 
^^^^ (d'iv + ip{v)dtv + f{v) = aw 
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in a bounded smooth domain $7 C M", n < 3. Here {v,w) = {v{t,x),w{t,x)) are unknown 
functions, A^^. is a Laplacian with respect to x, a > is a given parameter and / and ip are given 
nonlinearities which are assumed to satisfy some natural assumptions formulated in Section [2j 
Our interest to that problem is motivated by the following system arising in the mathematical 
biology: 

{dtu = j35w — ^{v)u — fu, 
dtv = fu- hv, 
dtw — dAxW + (3w = av, dnw\Q^ = 0, 

where a, /3, 6, d, /, h are given positive parameters and ^{v) is a given nonlinearity. This system 
has been introduced in [I4j in order to describe the development of a forest ecosystem (the 
unknown functions u, v and w are the densities of yang trees, old trees and the seeds respectively 
and the given nonlinearity 'y{v) describes the mortality of yang trees in dependence of the 
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density of the old ones) and has been recently studied analytically and numerically in [6l [3 [16] . 
Expressing u = {dtv + hv)/ f from the second equation of (jl.2p and insert it into the first one, 
we end up with the system of the type (|l.ip with respect to the variables (v^w). 

The main aim of the present paper is to study the long-time behavior of solutions of (jl.ip 
using the ideas and methods of the attractors theory. From the mathematical point of view, 
the problem considered is a coupled system of a second order ODE with a linear PDE (heat- 
like equation). Heuristically, it is clear that the dynamics of such coupled dissipative systems 
should depend drastically on the monotonicity properties of the ODE component. In the case 
where this ODE is "monotone", i.e., it cannot produce the internal instability (and all of the 
instability is driven by the coupling with the PDE component), one expects the asymptotic 
compactness and the existence of a smooth finite-dimensional global attractor with "good" 
properties. In contrast to that, in the non-monotone case, the ODE-instability may produce the 
asymptotic discontinuities and even may completely destroy the initially smooth spatial profile. 
Thus, in that case, the smoothing effect from the PDE component is not strong enough in 
order to suppress the development of discontinuities provided by the internal instabilities of the 
ODE-component and, as a result, spatial discontinuities and extremely complicated (in a sense, 
pathological) spatial structures may appear (see, for instance, [l] and [H] for the analysis of 
similar effects in the ID finite visco-elasticity) . We also mention that, although the existence and 
uniqueness of a solution of (II. 2p has been rigorously proved in the above mentioned papers [BJ [7] , 
very few has been done concerning the asymptotic behavior of solutions as t ^ oo. To be more 
precise, different types of w-limit sets of a single trajectory were considered there (associated 
with the different choice of the topology in the phase space) and their simplest properties were 
formulated, but the even the question whether or not they are empty remained open. As we will 
see below, some of them are indeed empty for the most part of the trajectories if the monotonicity 
assumption is essentially violated, see Remark 14.51 

In the present paper, we justify the above heuristics in a mathematically rigorous way on the 
example of the ODE-PDE coupled system (jl.ip . In particular, we show (in Section [3|) that the 
monotonicity arguments work perfectly if 

(1.3) f{v) >Ko>0, V £R. 

In this case, problem (jl.ip possesses indeed a smooth global attractor A of finite fractal dimen- 
sion in the proper phase space $oo- Moreover, due to the presence of a global Lyapunov function, 
this attractor generically has very nice properties (it is the so-called regular attractor in the ter- 
minology of Babin and Vishik [3j). Namely, it is a finite collection of the finite-dimensional 
unstable manifolds associated with the equilibria: 

(where TZ denotes a (generically finite) set of equilibria of problem (jl.ip and is an unstable 
manifold associated with the equilibrium uq € TZ, see Section [3]). Moreover, every trajectory of 
(jl.ip converges exponentially to one of that equilibria. We also mention that the first equation 
of (jl.ip is a second order ODE and, therefore, the monotonicity of / does not automatically 
imply the absence of the internal instability. For instance, the ODE 

y" + ^(y)y' + f{y) = h{t) 

may produce the non-trivial dynamics even if / is monotone and ip is strictly positive, say for 
the case of a given time-periodic external force h. By this reason, our proof of the monotonicity 
of the ODE component is based on rather delicate arguments related with the existence of the 
global Lyapunov function and associated dissipative integrals, see Section [3l 

The case where the monotonicity assumption is violated (which is considered in Section H] 
occurs to be (as predicted by the heuristics) much more complicated. In contrast to the monotone 
case, there is a very few hope to develop a reasonable global attractor theory here (no matter 
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in a strong or weak topology of the phase space), since, as a rule, even the equilibria set TZ 
is already not compact in the strong topology of the phase space and not closed in the weak 
topology. In addition, we see indeed a huge (uncountable) number of well-separated essentially 
discontinuous equilibria here, see Section HI 

Nevertheless, in the particular case of small coupling constant a, we succeed to give a complete 
description of the equilibria set TZ = TZa and verify that every trajectory of (jl.ip converges as 
t ^ oo to one of that equilibria. We mention that the standard Lojasiewicz technique for proving 
the stabilization seems to be non-applicable here even in the case of analytic non-linearities, since 
the equilibria set is not compact in any reasonable topology and the alternative technique of 
m and [19] (see also |13] where the pointwise stabilization for the non-smooth temperature 
driven phase separation model is proved) also does not work here since it is essentially based on 
the fact that the corresponding non-monotone ODE is a first order scalar ODE and cannot be 
generalized to the case of higher order equations. 

By this reason, we develop a new method of proving the stabilization, based on the theory of 
non-autonomous perturbations of regular attractors, see Appendix. 

Mention also that, as pointed out in the solutions with discontinuous densities are rather 
expected in a view of the forest ecosystem and a curve in Q where the density has discontinuities 
is called ecotone boundary. However, as our result shows, this "curve" is typically not smooth 
(and even not continuous) and may have an extremely complicated structure. 

The paper is organized as follows. Section[2]is devoted to the study of the analytical properties 
of problem (jl.ip such as existence and uniqueness, dissipative estimates in different norm, etc. 
The case of monotone nonlinearity / is considered in Section [3l in particular, the existence 
of smooth regular attractor is proved here. In Section H] we deal with the non-monotone case 
and, in particular, prove here the above mentioned stabilization result for the weakly coupled 
case. Finally, the Appendix is devoted to the derivation of the key estimate for our stabilization 
method which is, in turns, based on the perturbation theory of regular attractors. 

To conclude, we mention that it would be interesting to consider the regularization of (jl.ip 
in a spirit of a damped wave equation with displacement depending damping (see [9l [T71 [18] ) : 



We consider the following coupled system of a second order ODE with a heat equation: 



in a bounded 3D domain C M'^ with a smooth boundary. Here, {v,w) = {v{t, x),w(t, x)) are 
unknown functions, A^; is a Laplacian with respect to the variable x, a > is a given constant 
and (f and / are given nonlinearities, which satisfy the following assumptions: 




2. A PRIORI ESTIMATES, EXISTENCE AND UNIQUENESS 



(2.1) 



d'^v + ip(v)dtv + f{v) = aw, v{0) = vq, dtv{0) = v'q 
dtw - A^w + w = v, dnw\g^ = 0, If I^^Q = Wo 




(2.2) 
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V{v) > /3o > 0, 
f{v)v > -C + jo\v 
f'iv) > -K 



2+5 



(2.5) -Civit),w{t)) = -2i^ivit))dtv{t),dtv{t)) - a\\dtw{t)\\l,. 
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The aim of that section is to estabhsh a number of basic a priori estimates for that system which 
will allow us to verify the existence and uniqueness of a solution and to study its behavior as 
t ^ oo. We start with the following lemma which gives the global Lyapunov function for that 
problem. 

Proposition 2.1. Let the above assumptions hold and let {v{t),w{t)) be a sufficiently regular 
solution of problem ()2.ip . Introduce afunctional 

(2.4) C{v,w) := \\dtv\\l2 + 2{F{v),l) - 2a{v,w) + a\\V,,wf + a\\w\\l2, 

where F{v) := Jq f{s)ds and (•,•) is used for the inner product in . Then, the following 
equality holds 

d 

lit' 

Indeed, multiplying the first and the second equations of (|2.ip by dtv and adtw respectively, 
taking a sum and integrating over 0, we arrive at (|2.5p . 

Corollary 2.2. Let the above assumptions hold and let {v{t),w{t)) be a solution of (j2.ip . Then, 
the following estimate holds: 

(2.6) \\dMt)\\L2 + \\v{t)\\L2 + \\w{t)\\Hi < Q{\\{v,dtv,w)UJ 
for some monotone function Q independent of t and the solution 

Indeed, according due to our dissipativity assumption ()2.2p .3. 

(2.7) F(n) > -Ci + 7i|n|2+'^ 

for some new constants Ci and 71. Using this inequality, we easily check that 

(2.8) 72(||5t^|li2 + \\v\\l2 + IklllO - C2 < C{v,w) < Q{\\{v,dtv,w)UJ 

for some constants 72, C2 > and some monotone function Q. Integrating now equation (12. 5p 
by t and using that ip(v) > and a > 0, we arrive at (|2.6p . 

The next corollary gives the L^-dissipation integral for that problem. 

Corollary 2.3. Let the above assumptions hold. Then, 



(2.9) / \\dtv{t)\\l2 + \\dtw{t)\\l2 dt < Q{\\{vo,v'„wo)UJ 
Jo 

for some monotone function Q. 

Indeed, this estimate is an immediate corollary of (12. Sp . (12. Sp and the assumption that 93(f) > 
/3o >0. 

We are now going to verify that the solution is globally bounded in $00. 
Proposition 2.4. Let the above assumptions hold. Then, the following estimate is valid: 

(2.10) \Ht)\\L^ + \\dtv{t)\\L^ + \\w{t)\\L^nm < Qi\\ivo,v'o,wo)UJ 
for some monotone function Q independent of t and the solution. 

Proof. We first establish the L°°-bound for the tt;-component. Indeed, according to Corollary l2.2l 
the right-hand side v of the second equation of (j2.ip is bounded in L°°(R+, L^(r2)). Consequently, 
the standard regularity result for the heat equation gives 

(2.11) |k(t)||L- <C\\wm\L^e-' + C\\v\\L^^^^^L2)<Q{\\{vo,v'o,wo)UJ. 

(here we have implicitly used the restriction on the space dimension). 

Thus, we only need to establish the L°°-bounds for the t;-component. To this end, we will 
use the L°°-bounds for the if-component obtained before and will consider the equation for 



the u-component as an ODE for every (almost every, being a pedant) fixed x G il. Indeed, let 
y(t) := v{t,XQ). Then, this function solves 

(2.12) y"{t) + ip{y)y' + f{y) = h{t) = h^,^,{t) ■.= aw{t,xo). 
Multiplying this equation by y' + ey, we have 

(2.13) [{y' f + 2F{y) + 2eyy' + 2eR{y)]' + 2{^{y) - e){y' f + 2ef{y)y = 2h{y' + ey), 

where R{y) := Jq (p{s)s ds. Using that ip{v) is strictly positive and / is dissipative, we deduce 
from this equation that, for sufficiently small e > 

(2.14) j^S{y, y') + 7((y')' + v') < C(|/i(i)l' + 1), 
where 

S{y, y') ■■= {y'f + 2F{y) + 2eyy' + 2eR{y). 
and 7 is positive. Moreover, 

£o(2/' + {y'?) -C< S{y, y') < Q(y2 + (y'f) 

for some positive eo and C and monotone Q. Applying the Gronwall lemma to inequality (j2.14p . 
we conclude that 

(2.15) y{t)' + [y'it))' < Q{y{0)^ + y'{0)^ + \\h\\l^^^J 

for some monotone function Q which is independent of t and y, see eg, [1]. Taking now the 
supremum with respect to all xq (z and using (12. lip for estimating h, we deduce estimate 
(|2.10p and finish the proof of the proposition. □ 

We are now ready to verify the existence and uniqueness of a solution for the problem ()2.ip . 

Definition 2.5. A pair of functions {v{t),w{t)) is a solution of problem (12. ip if 

{V{t),dtv{t),w{t)) e $oo 

for every t >0 and (|2.ip is satisfied in the sense of distributions. 

Note that, from the first equation of (j2.ip . we see that dfv{t) € L°°{Q). Therefore v{t) G 
14^^'°°([0, r], L°°(il)) and the initial data for v is well-defined. Analogously, the w-component is 
continuous as a function with values, say, in Lp'iVi) and the initial data is again well-defined. 

Theorem 2.6. Let the above assumptions hold. Then, for every [vq,v'q,wq) G ^oo, problem 
()2.ip possesses a unique solution in the sense of Definition \2.5\ and this solution satisfies esti- 
mate (j2.10p . Moreover, any two solutions {vi{t),wi{t)) and {v2{t),W2{t)) satisfy the following 
estimate: 

(2.16) \\{vi{t),dmit),Wi{t)) - {v2{t),dtV2{t),W2{t))Uoo < 

< Ce'''\\{vi{0),dtvii0),wi{0)) - {v2{0),dtV2{0),W2mU^, 
where positive constants C and K depend only on the norms of the initial data. 
Proof. Let us first verify the uniqueness and estimate (j2.16p . Indeed, let 

{V{t),w{t)) := {vi{t),Wl{t)) - {v2{t),W2{t)). 

Then, these functions solve 

id'^v + (p{vi)dtv + [f{vi) - (p{v2)]dtV2 + [f{vi) - f{v2)] = aw, 
I dtw — + w = v. 



(2.17) 
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Multiplying now the first equation of (j2.17|) by dtv + ev, e > is a small positive number, using 
the fact that Vi,dtVi are globally bounded in L°° and applying the Gronwall inequality in a 
standard way (without integration by x\), we conclude that 

(2.18) \\dtv{t)\\l^f^^^ + \\v{t)\\l^ < 

<Ce'<\\\dMO)\\l^^n) + \Hml^in)) + cj\''^'-'^M^^^ 

for some positive constants C and K depending only on the L°°-norms of Vi and dfVi. Further- 
more, due to the maximum principle for the heat equation, we have the estimate 

(2.19) ||tf^(t)||L- < e"*||z/;(0)||Loo + f e~(*"') ||?;(s)||loo ds. 

Jo 

Inserting this estimate into the right-hand side of (|2.18p . we arrive at 

(2.20) \\dMt)\\l^^^^ + \\v{t)\\l^ < 

< C'e^\\\dMO)\\l^in) + imWU^n) + 11^(0) llioo ) + C e^(*-^) |lKs)||i^(^) ds. 
Applying again the Gronwall inequality to that relation, we conclude that 

(2.21) wdMmu^n) + Hmio^ < c'e'^widtvml^in) + hmlo^in) + Mml^)- 

This estimate, together with ()2.19p . give the desired L°°-estimate for the triple {v,dtv,w). In 
order to finish the proof of estimate ()2.16p . it remains to note that the desired estimate H^- 
norm of the w-component is immediate, since the L°°-control for the right-hand side of the heat 
equation for w is already obtained. Thus, the uniqueness and Lipschitz continuity ()2.16p are 
proved. 

So, we only need to prove the existence of a solution. It can be done in a standard way, based 
on a priori estimate ()2.10p . using the Banach fixed point theorem for proving the existence of 
a local solution and estimate (j2.10p for extending this solution globally in time, see eg, [12] for 
the details. □ 

Our next aim is to establish the basic dissipative estimate in the phase space •I'oo- 

Theorem 2.7. Let the above assumptions hold. Then, a solution {v{t),w{t)) of problem (12. ip 
satisfies the following dissipative estimate: 

(2.22) \Mt),v'{t),w{t))U^ <Q{\\{vo,v'„wo)UJe-^' + C, 
for some positive constants (3 and and monotone function Q. 

Proof. As we see from the proof of the previous proposition, the only problem is to obtain a 
dissipative estimate for the L^-norm of v{t). Indeed, if this estimate is obtained, analyzing the 
equation for the w-component analogously to (12.11j) . we deduce the dissipative estimate for the 
L°°-norm of w{t). This, in turns, gives the dissipative estimate for the right-hand side h{t) of 
(I2.12P and the Gronwall lemma applied to inequality (12.14p will finish the derivation of estimate 

So, we only need to obtain the dissipative estimate for the L^-norm. To this end, we multiply 
the first equation of (j2.ip by 2{dtv + ev), e > is a small number, which will be fixed below, 
and integrate over x € 0, after that we multiply the second equation of (j2.ip by 2a{dtw + 
£w), integrate over x E and take a sum of these two equations. Then, after the standard 
transformations, we end up with 

(2.23) ^Z{t)+2a\\dtwf+2{{^{v)-e)dtv,dtv)+2£f{v).v+2ae{\\VM\l2 + \Ml2) =4:aeiv,w), 



where 

(2.24) Z{t) := \\dtv{t)\\l, + 2{F{v{t)), 1) + 2s{R{v{t)), 1) + 2e{v{t),dtv{t))- 

- 2a{v{t),w{t)) + aiWVMmh + {I + e)\\w{t)\\l,). 

We now fix e > so small that 

emv{t)),l)\<l 

(it is possible to do due to estimate (j2.10p . of course, e will depend on the norm of the initial 
data). Then, due to (12. 71). we have 



(2.25) f32[\\dtv\\l2 + + \Mh + 1)] -C2< Z{t) < 

< /5i[l|5t^;|li2 + W^Mh + Mh + i\Hv)l 1)] + Ci, 

where the positive constants Ci and /3i are independent of e — )• and {v,w). Moreover, due to 
the fourth assumption of (|2.2p . 

(2.26) F{v) < f{v).v + Kv'^/2. 

Inserting estimates ()2.25p and ()2.26p into (|2.23p and using again the third assumption of ()2.2p . 
we deduce the differential inequality: 

(2.27) dtZ{t) + f3eZ{t) < Ce, 

where e depends on the norm of the initial data, but the positive constants (3 and C are inde- 
pendent of V and w. Integrating this inequality, we arrive at 

(2.28) Z(t)<[Z(0)-^]e-^^* + |. 

We see that, although the rate of convergence to the absorbing ball depends on the initial 
data (through the choice of e > 0), the radius of the absorbing ball is independent of e and, 
consequently, is independent of the norm of the initial data. This observation, together with 
estimate (j2.25p implies that 

(2.29) \\dtv{t)\\L2 + \\v{t)\\L2 + \\w{t)\\H^ < Qi\\ivo,v'o,wo)UJe-^' + C, 

for some positive 7 and C and a monotone function Q which are independent of t, f and w. 
Thus, the desired dissipative estimate in is obtained and Theorem 12.71 is proved. □ 

We now formulate several auxiliary results on the smoothing property for the tt;-component 
and the existence of dissipative integrals in stronger norms which will be essentially used in the 
next sections. 

Proposition 2.8. Let the assumptions of Theorem (|2.6p hold. Then, 'w{t) G T^^'^(r2) and 
dtw{t) G W'^'^{i^) for any t > and any p < 00 and the following estimate is valid: 

(2.30) \\w{t)\\w2.P^n) + \\dtw{t)\\w2,P^n) < + t''')Qp{\\ivo,v'o,wo)UJ 

for some positive exponent N and some monotone function Qp (depending only on p). 
Proof. Indeed, due to the smoothing property of the heat equation, the solution 9{t) of 

(2.31) dt9-A,9 + 9 = 0, el^^ = wo 
satisfies the following estimate 

(2.32) \\9{t)\\w^-pin) + \\dte{t)\\w^,,(^n) < Cpt-^\\wo\\m 

for some exponent N and positive constant Cp depending only on p, see eg, [15]. The remainder 
z{t) := w{t) — 9(t) solve the heat equation with zero initial data 

dtz - AxZ + z = v{t), ^Ij^Q = 0. 
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Moreover, using estimate (j2.10p together with the first equation of (|2.ip . we conclude that 
(2.33) HmL^in) + \\dtvmL^(^n) + l|5^(i)llL-{Q) < Q{\\{vo,v'o,wo)hJ. 

Using that estimate together with the PF^'^-regularity estimate for the heat equation, we arrive 
at 

\\^it)\\w^'P{n) + \\dtw{t)\\w2,p^^) < QpdK^o, ''^0,^^0)11*00) 
which, together with estimate (j2.32p . finishes the proof of the proposition. □ 

Proposition 2.9. Let the assumptions of Theorem \2.6\ hold. Then, the following stronger 
version of dissipative integrals exist: 

fl) /,°° \\dtw{t)\\l^^^ dt < Q{\\{voA,wo)hJ, 

[2) \dMt,xo)\^ + \dtv{t,xo)\^ dt < eT + CeQ{\\{vo,viwo)UJ, 

where e > is arbitrary, xq (z Q is almost arbitrary, Cg > depends only on e and Q is some 
monotone function. 

Proof. We first note that, due to Lemma 12.81 we may assume without loss of generality that 
dtw{0) € H^. Differentiating the equation for the w-component by t and denoting z := dtw, we 
get 

(2.35) dtz - A^z + z = dtv, z\^^^ = dtv{0). 
Applying the L'^-regularity theorem for that heat equation, we will have 

(2.36) r l|A.z(s)||^2 ds < C\\z{0)\\l, +C T \\dtv{s)\\l2 ds, 
Jo Jo 



(2.34) 



where the constant C is independent of T. Together with estimate (12. 9p and embedding C 
, it gives the desired first estimate of (j2.34p . Let us now prove the second estimate of (j2.34p . 
To this end, we multiply equation (I2.12P by 2y' (without integration by x!). This gives 

(2.37) {{y'f + 2F{y) - 2awy)' + 2^{y)(:y'f = -2aydtw{t, xq). 
Integrating this equality over t E [0,T], estimating 

\2aydtw{t,XQ)\ < C\\v{t)\\L^\\dtw{t)\\L^ < e + e^'^a-^\\v{t)\\loo\\dtw{t)\\loo , 

and using the first estimate of (I2.34p together with the strict positivity of ip and the fact that 
the L°°-norm of v is under the control, we deduce that 

(2.38) [ [y'{t)fdt<eT + e-^Qi\\vo,v'o,woUJ 

Jo 

for some (new) monotone function Q which is independent of T. This gives the second estimate 
of (I2.34P for the term dtv{t,xo). Thus, in order to finish the proof of the proposition, we only 
need to estimate the term dfv{t,xo). To this end, we differentiate the first equation of (12. ip by 
t and denote q{t) := dtv(t,xo). Then, we get 

q" + vivW + ^'{v)q^ + f'{y)q = adtw{t, xo). 

Multiplying this equation by 2q' , integrating by time and using that the L°°-norms of v, dtv and 
dfv are under the control, we arrive at 

'^[q'{t)fdt < Q{\\voA^w4^J [l + j\y\t)? + \\dtw(t)\\l^ dt 

for some monotone function Q which is independent of T. Inserting estimate ()2.38p and the first 
estimate of ()2.34p to that inequality (and scaling the parameter e if necessary), we obtain the 
desired control for the integral of d^v and finish the proof of the proposition. □ 
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We conclude this section by showing that, if the initial data {vo,v'q,W()) is smooth, the solution 
{v{t),w{t)) remains smooth for all t. 

Proposition 2.10. Let the assumptions of Theorem ()2.6p hold. Assume, in addition, that 
(2.39) {vo,v'o,wo)£W^'"^{n). 

Then, the solution {v{t),dtv{t),w{t)) G W^'°°{n) for any t > and the following estimate holds: 



Kt 



(2.40) ||v(t)|!t^i,oo + \\dtv{t)\\yi^i,oo + ||?i;(t)||^yi,oo < C\\{vo, v'q, wo)\\[wi'°^]3 
for some positive constants C and K (which depend on the -norms of the initial data). 

Proof. The desired estimate for the tt;-component is factually obtained in Proposition 12.81 thus, 
we only need to estimate the v component. Let xi and X2 be two arbitrary points of and let 
z{t) := v{t,xi) — v{t,X2). Then, this function satisfies the following ODE 

(2.41) Z"{t) + ipiv{t, Xl))z'{t) + [ip{v{t, Xl)) - ip{v{t, X2))]dtv{t, X2) + 

+ [f{v{t,xi)) - f{v{t,X2))] = /ixi,x2(i) := a{w{t,xi) - w{t,X2)). 
Multiplying this equation by dtz{t) and arguing exactly as in (I2.18p . we arrive at 

(2.42) \z{t)\^ + \z'{t)\' < C{\z{0)\^ + |z'(0)|2)e^* + C [\'"^'-'^\h^,,,,{s)\Us, 

Jo 

where the constants C and K depends on the L°°-norm of the solution. Furthermore, since the 
VF^'°°-estimate for the iM-component is already obtained, we have 

sup -. ^ f e''^'-'\h^,,^,{s)\^ds<Ci re^(*-^)||u;(s)||^i,oods<Ce^*. 

Dividing finally inequality (j2.42p by \xi — and taking the supremum over xi,X2 G from 
the both parts of the obtained inequality, we obtain the desired estimate for the Ty^'°°-norms 
of V and dtv and finish the proof of the proposition. □ 

Remark 2.11. Arguing analogously, one can show that if the initial data belong to C*^, the 
solution will be of class for every t > 0. Thus, the blow up in finite time of the higher norms 
cannot occur. However, there is a principal difference between estimates (|2.22p for the L°°-norm 
and estimate (|2.40p for the T^^'°°-norm of the solution. Indeed, the first estimate is dissipative 
and shows that the L°°-norm of the solution cannot grow and even gives the absorbing ball 
in that norm. In contrast to that, the VF^'°°-norm, a priori, may grow exponentially and, in 
this sense, the solution may become "less and less regular" as t — > 00 (i.e., it may tend to a 
discontinuous limit). As we will see in the next sections, the answer on the question whether or 
not it really happens depends in a crucial way on the monotonicity of the nonlinearity /. 

3. The monotone case: asymptotic compactness and regular attractor 

According to the results of the previous section, equation (12. ip is uniquely solvable in the 
phase space ^>oo and the solution operators 

(3.1) S{t){vo,v'Q,wo) := {v{t),dtv{t)Mt)) 

generate a dissipative semigroup in ^oa- The aim of this section is to study the long-time 
behavior of solutions as i — >• cx) in the particular case where the nonlinearity / is strictly 
monotone: 

(3.2) f'{v) > Ko > 0. 

As we will see, in that case, the associated semigroup is asymptotically compact and possesses 
a smooth global attractor A in ^oo- Moreover, due to the Lyapunov functional, this attractor 
can be described as a finite union of finite-dimensional unstable manifolds. Our proof of the 
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asymptotic compactness is based on the following lemma which can be considered as a refinement 
of estimate ([230]) • 

Lemma 3.1. Let the assumptions of Theorem \2. 1\ hold and let, in addition, ()3.2p he satisfied. 
Let us also introduce, for any h > the following (semi)norm on the space 

(3.3) ||f ||^;^^i,cx) := sup i i 

(being pedants, we would write esssup instead of sup). Then, every solution {v{t),w{t)) of 
problem (j2.ip satisfies 

(3.4) \\vm^^,^ + \\dtvit)\\^..o.<^e'^' + C2, 

where the positive constants (3, Ci depend on the ^oo-norms of the initial data, but are indepen- 
dent of t and /i — > 0. 

Proof. Analogously to the proof of Proposition 12.10] we introduce a function z{t) := v{t,xi) — 
v{t,X2) which solves equation ()2.4ip . But, using the monotonicity assumption ()3.2p and the 
dissipation integrals ()2.34p . we are now able to suppress the exponential divergence in estimate 
()2.42p . To this end, we multiply equation p.4ip by z' and transform the term containing the 
nonlinear ity / as follows: 

(3.5) [f{v{t, xi)) - f{v{t, X2W = l/2[/(t)z2]' _ 1/2V {t)z\ 

where l{t) := f'{sv{t, xi) + (1 — s)v{t, X2)) ds > kq > ^ and its derivative can be estimated as 
follows: 

(3.6) \l'{t)\ < C{\dtv{t,xi)\ + \dtv{t,X2)\), 

where the constant C depends on the L°°-norm of the initial data, but is independent of t and 
Xi. Then, using the positivity of (p and the L°°-bounds for v, we get 

l/2((z')2 + /^7 + 7(^')'<C|/i.„.,P + C7(|ait;(t,xi)| + |ait;(t,X2)|)z2 

for some positive constants 7 and C. Multiplying now equation ()2.4ip by ez (where e > is a 
sufficiently small positive number) and taking a sum with the above inequality, we infer 

(3.7) l/2((z')^ + Iz^ + 2ezz'y + (7 - e){z' f + elz^ < 

< c\h^u^,\^ + ci\dtv{t,xi)\ + \dtv{t,x2mz^ + (z'f)- 

Let now 

C^{t) := {z'f + lz^ + 2£zz' . 
Then, since l{t) > kq > 0, we may fix e > to be small enough that 

(3.8) k(z2 + [z'f) < C,{t) < Ki{z^ + (z'f) 

for some positive k and ki. This, inequality, together with the evident estimate |x| < /3 + /3~^x'^ 
allows to transform (13.71) to 



(3.9) ^A(t) + (7 - C{\dtv{t,xi)\^ + \dtv{t,X2)f))C,{t) < C\h.„xM^- 

Applying the Gronwall inequality to this relation, we arrive at 

C,{t)<C,{0)e-^o^^^)<i^ + C r e-I^''^'^'''\h,,,,,{t)\Ut 

Jo 
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2 , \a ^\i2\ 



i: 



withi^(s) := ^—C{dtv{t,xi)\ +\dtv(t,X2)\ )■ Using the dissipation integral (|2.34p with e := 7/2, 
we see that 

r-T 

K{s)ds>j{T-t)/2-C, 

It 

where the constant C depend on the $oo-iiorm of the initial data. This estimate, together with 
the bounds (I3.8p . give the non-divergent analogue of estimate (12.42P : 

(3.10) + \z'{t)\^ < C(|z(0)p + |z'(0)|2)e-^*/2 + C f e-<'^'^'^'^\K,,.,{s)\^ ds 



Jo 

for some positive C and 7 depending only on the <l>oo-iiorm of the solution. 

In order to deduce the desired estimate (|3.3p from (|3.10p . we note that, due to Proposition l2.8l 
we may assume without loss of generality that ||?i;(t)||54/i,oo < C for alH > and, consequently, 

sup -. ^ f e-^^'~'y^\h,,,^,{s)f ds < a f e-^'^'-'^/^\\w{s)\\l,^,^ ds < 

\xi-X2\>h — ^2] Jo Jo 

Moreover, obviously. 

Dividing now inequality (|3.10p by |xi — and taking the supremum over all Xi G Q, |xi — X2I > 
h, we deduce the desired estimate (j3.3p . Lemma [STT] is proved. □ 



Our next step is to verify the existence of a global attractor A for semigroup ()3.ip associated 
with problem (12. ip . We recall that, by definition, the global attractor A should satisfy the 
following properties: 

1) A is compact in <I'oo; 

2) A is strictly invariant: S{t)A = A; 

3) A attracts the images of all bounded sets as t — )• 00, i.e., for any bounded set B C $00 and 
any neighbourhood 0{A) of A in ^qo, there exists time T = T{B,A) such that 

S{t)B C 0{A), for t > T. 

We also recall that the attraction property can be also reformulated in terms of the non- 
symmetric Hausdorff distance between sets in ^oo- 

(3.11) lim disiiS{t)B,A) = 0, 



see eg, p] for the details. 

Theorem 3.2. Let the assumptions of Lemma \3.1\ hold. Then, the semigroup S{t) associated 
with problem ()2.ip possesses a global attractor A in the phase space $00 • This attractor is 
bounded in [W^''^{Q,)]^ and has the following structure: 



(3.12) a = k:\ 



t=0' 



where IC C C{,(M, <I>oo) is a set of all solutions of problem (12. ip which are defined for all t M 
and are globally bounded. 

Proof. In order to deduce the existence of a global attractor from Lemma 13. H we will use the 
so-called Kuratowski measure of non-compactness. Recall that, by definition, the Kuratowski 
measure of non-compactness a{B) of a set B is infinum of all r > for which it can be covered 
by a finite number of r-balls, see e.g., JT] for details. To be more precise, we need the following 
lemma. 



12 



M.EFENDIEV AND S.ZELIK 



Lemma 3.3. Let 

B ■.= {v e L°°(0), II^IIloo + ||t;||^^i.oo < R} 

for some R and h > 0. Then, its Kuratowski measure of non- compactness of the set B can be 
estimated as follows: 

(3.13) a{B) < Rh. 

Proof. Let Sh be the standard averaging operator 



{Shv){x) := / Dhix,y)v{y)ds, 

where the smooth non-negative kernels Dh{x,z) are such that 

1. supp L'/j(x, •) C {z S il, \z — x\ < h} 



(3.14) 



2. J^^Dh{x,y)dy = 1, 

2,. \Dh{x,y)\ + \VxDh{x,y)\ <Ch, x,y e 



(since is assumed to be smooth, such kernels exist). 

Let also B^ := Sh{B). Then, on the one hand, the set Bfi consists of smooth functions and, 
in particular, is bounded in C^{Q). By the Arzela-Ascoli theorem, it means that Bh is compact 
in L°°{n). 

On the other hand, 

\{Shv){x) - v{x)\ < / Dh{x,y)\v{y) - v{x)\dy < \\v\\y^^l,^h Dh{x,y) dy < Rh. 



h 



Thus, B C Bfi + Rh and B^ is compact. This gives estimate (|3.13p and finishes the proof of the 
lemma. □ 

We are now ready to finish the proof of the theorem. Lideed, due to Proposition 12.81 we know 
that the w-component is bounded in W'^'^{Q,) for every t > and, consequently, the w-component 
of S(t)B is precompact in L°° n for any bounded set B. So, the Kuratowski measure of 
non-compactness for S{t)B is determined by the u-component only. Moreover, Lemma 13.11 
guarantees, that 

v{t),dtv{t) C {n G L°°(0), ||n||Loo + ll-uH^^i.cx, < R} 
ii t > T{h) is large enough (but R is independent of h). This, gives that 
(3.15) lim a{S{t)B) = 

for any bounded set B. 

Since the semigroup S{t) is Lipschitz continuous with respect to the initial data (see Theorem 
12. 6p and dissipative (see Theorem 12. 7p , the convergence of the Kuratowski measure (I3.15j) to 
zero implies the asymptotic compactness of the semigroup S{t) and the existence of a global 
attractor A, see [llj. The structure (j3.12p of the attractor is also a corollary of that abstract 
theorem and the fact that A is bounded in W^'°° follows from estimate (j3.3p (together with the 
fact that the constant C2 is independent of h). Thus, Theorem 13.21 is proved. □ 

Our next task is to establish the regular structure of the attractor A provided by the Lyapunov 
functional. To this end, we need to make some preparations. As a first step, we establish the 
differentiability of the semigroup S{t) with respect to the initial data. 

Proposition 3.4. Let the assumptions of Theorem \2.6\ hold. Then, the associated semigroup 
S{t) is Frechet dijferentiable with respect to the initial data for every fixed t and it's Frechet 
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derivative D^[S{t)$] G C{^oo, ^oo) is Lipschitz continuous with respect to the initial data ^ € <i>oo 
and the following estimate holds for every bounded set B C $00 • 

(3.16) ||5(t)||ci,i(B,$^)<Ce2^*, 
where the constants C and K depend only on the radius of B. 

The proof of this proposition is straightforward and standard, so, in order to avoid the tech- 
nicahties, we rest it to the reader. 

At the next step, we need to study the equilibria of problem (I2.ip . 

Proposition 3.5. Let the assumptions of Theorem \3.^ hold. Then, any equilibrium {vq,Wo) € TZ 
(the set of all equilibria) of problem (j2.ip solves the following semilinear elliptic equation: 

(3.17) - A^w + w = f-\aw), 9„|^^ = 0, v = {-A^ + l)w 

(f~^ exists since f is now assumed to be monotone). Moreover, the equilibrium {vq,wo) is 
hyperbolic if and only if wq is hyperbolic as a solution of ()3.17p . i.e., if the equation 

(3.18) -AJ + e= [f-^]' {awo)a9 

has only trivial solution = 0. In particular, for generic f , all of the equilibria {vq, wq) € TZ are 
hyperbolic and TZ is finite. 

Proof. Indeed, the equations on equilibria for problem (j2.ip 

f{v) = aw, —AxW + w = V 

are equivalent to (j3.17p . Let us verify the assertion on hyperbolicity. Indeed, the asymptotic 
compactness of the semigroup S{t) implies in a standard way that the essential spectrum of the 
operator D^S{1) lies strictly inside of the unit circle. Thus, only eigenvalues of finite multiplicity 
are possible on the unit circle. Any such eigenvalue generates a time-periodic solution {z, 6) of 
the associated equation of variations 

(dfz + ip{vo)dtz + f{vo)z = a9, 
I dtO - A J + 9 = z. 



(3.19) 



However, analogously to the nonlinear problem (j2.ip . the linearized problem (j3.19p possesses a 
global Lyapunov function (in order to find it, one needs to multiply the first and the second 
equations by dtz and adtO respectively, take a sum and integrate over O). Thus, every periodic 
solution of that linearized problem must be an equilibrium: 

f'{vQ)z = ad, z = -AJ + e 

and, consequently, z must solve (|3.18p . Vise versa, any nontrivial solution z of (j3.18p generates 
a non-trivial equilibrium of (j3.19p by setting z = —A^O + 9. 

Finally, the last assertion that generically TZ is finite and all of the equilibria are hyperbolic 
is a standard corollary of the Sard theorem, see eg, [3]. Proposition 13.51 is proved. □ 

Thus, we will assume from now on that all of the equilibria [vq, wq) (zTZ are hyperbolic (which 
automatically implies that TZ is finite). Furthermore, for any ^0 '■= {vo,wo) G TZ, we define the 
associated unstable set by the usual expression 

(3.20) Mt ■■= {{vo,Vq,wo) G $00, 3{v(t),w{t)) which solves t^J^i for t < such that 



{v{0),dtv{0),w{0)) = {vo,Vq,wo) and lim {v{t),dtv{t),w{t)) = {vo,Vq,wq)}. 



In other words, the unstable set consists of all complete trajectories of (j2.ip which stabilize 
to (,0 as t ^ —00. 
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It is well known (see e.g., [3]) that, for hyperbolic equilibrium ^ the set is locally 
(near ^o) a finite-dimensional submanifold of $00 and its dimension equals to the instability 
index of ^o- But, in order to prove that the whole is a submanifold of $00 one needs 
the semigroup S{t) to be injective (in other words, problem (|2.ip should possess the so-called 
backward uniqueness property, see again [3]). 

Proposition 3.6. Let the assumptions of Theorem \2. 61 hold. Then, the semigroup S{t) associ- 
ated with equation (j2.ip is injective, i.e., the equality S{T)(,i = 5'(T)^2; for some T > 0, implies 
that ^1 = ^2 • 

Proof. Indeed, let {vi{t),wi{t)) and {v2{t),W2{t)) be two solutions of problem (j2.ip and let 
(z{t),w{t)) be their difference. Then, these functions solve (j2.17p . Let us rewrite this equation 
in the form 

(3.21) dt^ + B^ = Vm, 

where ^(t) := iz{t),dtz{t),e{t)), 

1 1 0' 

B={0 1 |, Vit) := {-[l^{t)dtV2{t) + lf{t)] l-ip{vi{t)) a 

1 0, 





and 



'1 i-i 
l^{t):= I ip\svi{t) + {l- s)v2{t))ds, lf{t):= if' {svi(t) + {I - s)v2{t)) ds. 



Let us consider equation (I3.2ip in a Hilbert space H := [L^(r2)]'^. Then, obviously, i3 is a positive 
self-adjoint (unbounded) operator in H and the operator 'P(t) is uniformly bounded for all t > 0. 
Then, the classical backward uniqueness theorem of Agmon and Nirenberg (see [2j) is applicable 
here and, consequently, ^(T) = implies that ^(0) = 0. Proposition 13.61 is proved. □ 

Corollary 3.7. Let the assumptions of Theorem \3.2\ hold. Then, for any hyperbolic equilibrium 
^0 := {vq,wo) of problem (j2.ip . the associated unstable manifold Ai'^^ is a finite- dimensional C^- 
submanifold of diffeomorphic to M^, where N is the instability index of the equilibrium ^o- 

This result is a standard corollary of the existence of a Lyapunov function, injectivity and 
smoothness of the semigroup S{t), see [3]. 

We are now ready to formulate a theorem on the regular structure of the attractor A which 
can be considered as the main result of the section. 

Theorem 3.8. Let the assumptions of Theorem \3.2\ hold and let, in addition, all of the equilibria 
S,o & TZ be hyperbolic. Then, 

1) any non- equilibrium complete trajectory i{t), t € M of the semigroup S{t) belonging to the 
attractor is a heteroclinic orbit between two different equilibria ^_ and 

(3.22) lim e(t) = C±, 

t±oo 

where ^_ 7^ 

2) The attractor A is a finite union of finite- dimensional submanifolds o/<5oo-' 

(3.23) A = U^,enMl, 

where Al^ is an unstable manifold of ^0 ^ 

3) The rate of attraction to A is exponential, i.e., for any bounded set B C 

(3.24) disi{S{t)B,A) < Q{\\BUJe-^' 
for some monotone function Q and positive constant 7. 
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Indeed, this theorem fohows from the abstract result of Babin and Vishik on regular attractors, 
see [3j (see also [20j) and Propositions 12.11 13.41 and 13. 6i 

Remark 3.9. Theorem 13.81 shows that the long-time behavior of solutions of problem (12. ip is 
"extremely good" if the monotonicity assumption ()3.2p holds. As we will see in the next section, 
this behavior is "extremely bad" if the monotonicity assumption is essentially violated. 

4. The non-monotone /: stabilization for the case of weak coupling 

The aim of this section is to understand how the asymptotic behavior of (j2.ip may look like 
when the monotonicity assumption (j3.2p is violated. To this end, we will consider below the 
case of small coupling constant a, where the dynamics is, in a sense, determined by the limit 
ODE 

(4.1) y"it)+^iy{t))y'{t)+f{yit)) = 0. 

In particular, in that limit case, the value of v{t, x) at x = xq is determined by the value 
of {vq{x),Vq{x)) at X = Xq Only and, consequently v{t,x) and v{t,y) evolve independently if 
X ^ y. Thus, if (j4.ip has more than one equilibrium, the most part of trajectories will tend to 
a discontinuous equilibria, no matter whether or not the initial data {vq,Vq) is continuous. As 
we will see, the same property preserves for the case of small positive coupling constant a. 

To be more precise, we assume that the limit equation (14. ip possesses a regular attractor in 
M^, i.e., that 

(4.2) /'(tio) / 0, for ah uq such that /(uq) = 

(since the existence of a global Lyapunov function and dissipativity are immediate, only the 
hyperbolicity assumption on the equilibria should be postulated). 

We start our exposition by verifying that the dissipative estimate (|2.22p is uniform with 
respect to a ^ 0. 



Proposition 4.1. Let the assumptions of Theorem 2.1 hold. Then the positive constants (3 and 
C^, and monotone function Q in the dissipative estimate (|2.22p are independent o/a — > 0. 

Proof. In order to verify this assertion, we need to check that the most part of estimates of 
Section [2] are uniform with respect to a — >■ 0. We start with estimate (12. 6p . From the first point 
of view (based on the form of the global Lyapunov function (12. 5p ). one may guess that it is 
non-uniform and only a|[?i;(t)||^2 is uniformly bounded. However, the Lyapunov function gives 
the uniform bound for the L'^-norm of the u-component. The standard L^-estimate for second 
equation 

(4.3) dtw - A^vu + UJ = v{t) 

of ()2.ip gives after that the uniform and L°°-bounds for the w-component. 

Thus, the uniformity as a —t' is verified for estimates (|2.6p and (12. lip . The uniformity of 
the L°°-bound (j2.10p follows from (j2.1ip exactly as in Proposition 12.41 

So, it only remains to check the dissipative estimate (I2.22P and, following the proof of Theorem 
2.71 we see that only the uniformity of the L'^-estimate (I2.29P is necessary. 

Analogously to (12. 6p the function Z(t) (defined by (I2.24p ) can be estimated as follows 

(4.4) MWdtvWh + «(l|V.u;||i2 + ||u;||i.) + (|F(z;)|, 1)] - C2 < Z{t) < 

< MWdtvWh + aiWVMlh + \Ml2) + mv)\, 1)] + Ci, 

where /3j and Cj are now uniform with respect to a — > 0. By this reason, estimates (j2.27p and 
(j2.28p do not give immediately the uniform analogue of ()2.29p . but only the uniform dissipative 
estimate for the L^-norm of v(t). Combining after that this estimate with the usual L^-estimate 
for the heat equation ()4.3p . we verify that estimate ()2.29p is indeed uniform as a ^ 0. Exactly 
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as in Theorem 12.71 this gives the uniformity of estimate (|2.22p and finishes the proof of the 
proposition. □ 

Thus, due to the previous proposition, the radius of the absorbing ball in $00 for problem 
()2.ip is uniform with respect to a — )• 0. In particular, the ^>oo-norm of any equilibria of that 
problem is uniformly bounded. Denoting the set of equilibria for problem (|2.ip by 7^^, we may 
conclude that 

(4.5) ||7^a||$^<C, 

where the constant C is independent of a. 

This observation together with the hyperbolicity assumption (I4.2p allow to give a complete 
description of the equilibria set TZa if a > is small enough. 

Proposition 4.2. Let the assumptions of Theorem \2.6\ hold and let, in addition, the limit 
hyperbolicity assumption (|4.2p be valid. Denote these hyperbolic equilibria by {ui,--- ,un}- 
Then, there exists > such that, for every a < and every partition 

(4.6) i7 = Oi U • • • Ul^AT 

on disjoint measurable sets: iljnilj = fori 7^ j, there exists a unique equilibrium {vo,0,Wo) G 
$00 of problem ()2.ip such that 

N 

(4.7) vo = vo + 9, vo{x) ■.= ^UiXn,{x), \\9\\l^ < Ca, 

1=1 

where xv{x) is a characteristic function of the set V and the constant C is independent of a. 
Moreover, every equilibrium {vq,0,Wo) € $00 can be presented in such form. 

Proof. Indeed, in order to find the equilibrium, we need to solve 

f{vo) = awo, A^wo -wo = vo 

which we rewrite in the form of a single equation on vq in L°°{^}): 

(4.8) /(^;o) = a(-A,. + ly^o- 
We note that the function 

F{v,a) := fiv)-ai-A^ + l)-^v 
belongs to C2(L°°(17) x ]R,L°°(0)). Moreover, its derivative 

D,F{do,0) := f'ivo) 

is invertible in L°° (due to the hyperbolicity assumption (14. 2p and the norm of the inverse 
operator is uniformly bounded with respect to the choice of a partition. In addition, F{vo, 0) = 0. 
Thus, the existence and uniqueness of the equilibrium vq in a small neighbourhood of vq if a is 
small follows from the implicit function theorem. 

Let us now verify that any equilibrium (vq,wo) can be presented in that form. Indeed, let 
00 be an arbitrary equilibrium. Then, according to (14. 5p . ||ifo||-L°° — C where C 
is independent of a. Therefore, 

\f{vo{x))\ <Ca, X E 

Since all of the roots f{z) = are hyperbolic, for sufficiently small a, we conclude from here 
that 

(4.9) vq{x) G Oca{uk{x)), X eVt 
for some root ti^j^^.) of f{z) = 0. Fixing now 

(4.10) := {x G Jl, k{x) = i}. 
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we see that the equUibrium {vo,wo) indeed has the form of (j4.7p and the proposition is proved. 

□ 

Remark 4.3. We see that, in contrast to the case of monotone /, we now have the uncountable 
number of different equihbria (ah of them are hyperbohc in $00) most of which are discontinuous 
(we have only finite number of continuous equihbria associated with trivial partitions of 0). 
Moreover, using the explicit description given in the previous proposition, it is not difficult to 
show that the set TZa is not compact in the strong topology of the space <I>oo and is not closed in 
the weak-star topology of that space. By this reason, the possibility to apply the strong/weak 
global attractor theory to that problem seems very problematic. However, as the next theorem 
shows, any trajectory {v{t),w{t)) still converges to one of the equilibrium from T^a 

as f — >■ 00. 

Theorem 4.4. Let the assumptions of Proposition \4-^ hold. Then, there exists ao > such 
that, for every a < ao every trajectory {v{t),dtv{t),w{t)) of problem (12. ip stabilizes as i —t- 00 
to some equilibrium {v,0,w) € TZa in the topology of U'{Q): 

(4.11) hm \\{v{t),dtv{t)Mt)) - {v,Q,w)\\^L,in)f = 

for any 1 < p < 00. 

Proof. The proof of that convergence is strongly based on the perturbation theory of regular 
attractors and Proposition [5T] (see Appendix). Indeed, due to Propositions 14. 1 1 and 12.81 we may 
assume without loss of generality that (i;(0), 9(w(0), t«(0) belongs to the absorbing ball Bn in 
<^oo with the radius R independent of a and that 

(4-12) ||w^||c6(M+xn) + \\dtw\\ci,{R+xn) < C, 

where the constant C is also independent of a. Thus, the first equation of (12. ip 

(4.13) dfv{t, x) + Lp{v{t, x))dtv{t, x) + f{v{t, x)) = aw{t, x) 

can be treated as an ODE for every fixed x G Moreover, due to the hyperbolicity assumption 
(j4.2p and uniform estimate (j4.12p the right-hand side of ()4.13p can be treated as small non- 
autonomous perturbation of the ODE 

(4.14) u" + ip{u)u' + f{u) = 0. 

Thus, the assumptions of Proposition 15.11 hold for problem ()4.13p for every fixed x G if a < ao 
for sufficiently small positive ao- Due to this Proposition, we have the estimate: 

(4.15) / \dMt,x)\ + \dtv{t,x)\dt <Ci + C2a [ \dtw{t,x)\dt, 



where the positive constants Ci and C2 are independent of T, a and x € fi. Integrating this 
inequality by x € ^2, we arrive at 

(4.16) r 119^(011^1 + \\dtv{t)\\Li dt < Ci\n\ + Caa / \\dtwit)\\L, dt. 
Jo Jo 

In order to estimate the integral into the right-hand side of (I4.16p . we differentiate the second 
equation of ()2.ip by t, denote 6 := dtw, multiply it by sgn6{t) and integrate over Q. Then, due 
to the Kato inequality, we arrive at 

(4.17) dtmt)\\Li + \\9\\L^ < \\dtv{t)\\Li. 
Integrating this inequality, we have 

\\dtw{t)\\L^ < \\dtwm\L^e-' + [\-'^'-'^\\dtv{s)\\Lids. 
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Integrating the obtained inequality once more over t € [0, T] and using tliat dtw{0) is uniformly 
bounded, we arrive at 

(4.18) r \\dtw{t)\\Lidt<C+ r \\dtv{t)\\Lidt, 

Jo Jo 

where C is again independent of a and T and the trajectory. 

Inserting M.lSh into the right-hand side of (I4.16p and assuming that a is small enough, we 

finally deduce the following L^-dissipation integral 

(4.19) r \\dMt)\\L^ + \\dtvit)hr + \\dtwit)h. dt < C, 

Jo 

where the constant C is independent of T. 

Thus, we have proved that {v{t),dtv{t),w{t)) converges to some ^ G $oo in the L-^(r2)-norm. 
Moreover, since we have the control of the L°°-norm, the interpolation inequality gives the 
convergence in for any p < oo: 

(4.20) lim \\iv{t),dtv{t),w{t)) - (\\[LPin)Y^ = 0. 

Thus, we only need to verify that ^ € TZa is an equilibrium. To this end, we will use the so- 
called trajectory approach (see [5] for the details) and consider positive semi-trajectories instead 
of points in the phase spaces. Indeed, arguing exactly as in the proof of estimate ()2.16p . but 
taking the L^-norm instead of the L°°-norm, we see that 

(4.21) \\{vi{t),dMt),Wi{t)) - {v2{t),dtV2(,t),W2m\lLPm^ ^ 

< Ce''\\{vi{0),dtvii0),wi{0)) - {v2iO),dtV2iO),W2m\\[LPin)]^. 
Define now the map § : $00 L°°(]R-(_, ^00) via the expression 

(4.22) S : (t;o,t;^,t/;o) ^ ivi-),dtvi-),wi-)) 

and let /C+ := §($00)- Then, estimate (|4.2ip (together with the obvious fact that JC+ C 
L°°(R+, <I>oo)) shows that the map S realizes a Lipschitz continuous homeomorphism between 
spaces $00 and /C+ endowed by the topology of [LP{Q)]^ and ^tr '■= L°°{R+, [LP{Q)]^) respec- 
tively. The solution semigroup S{t) is conjugated via that homeomorphism to the semigroup 
T{t) of temporal shifts on ^tr- 

(4.23) s(t) = or(t) oS, (r(t)O(s) :=e(t + s), e e /C+, t,s>0. 

Thus, the convergence ()4.20p implies that 

T{t)§^o :=T{t)§{v{0),dtv{0),w{0)) 

in the space ^tr- 

Let us first check that §^ € /C+, i.e., that the limit trajectory ^(t) solves equation (12. ip . In 
other words, we need to show that /C+ is closed in ^tr- 

To this end, we need to show that it is possible to pass to the ^t^-liiiiit in equations (|2.ip in 
the sense of distributions for any sequence (,n{t) := {vn(t) , dtVnit) , Wnit)) converging in ^tr to 
some ^{t) := {v{t),dtv{t),'w{t)) and bounded in L°°(R+,$oo)- Indeed, the passage to the limit 
in all linear terms are evident and only the passage to the limit in the nonlinear terms (p{v) and 
f{v) may a priori be problematic. But it is not the case, since convergence in ^tr implies the 
convergence almost everywhere (up to extracting a subsequence) and this allows us to conclude 
in a standard way that f{vn) f{v) and ip{vn)dtVn ip{v)dtv (here we have implicitly used 
that Vn is uniformly bounded in L°°). Thus, the limit function ^(t) solves indeed problem (|2.ip . 

We are now ready to verify that ^(t) is an equilibrium which will finish the proof of the 
theorem. Indeed, due to the dissipation integral (j2.9p . we see that 

\\dtv{s + •)IIl2{r+xQ) + \\dtv{s + ■)\\L^R+xn) + \\dtw{s + ■)\\L^R+xn) ^ 
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as s — )• (X). Thus, for the hmit functfon {v{t),dtv{t),w{t)), we have dtv = d^v = dtw = and ^ 
is indeed an equilibrium. Theorem 14.41 is proved. □ 

Remark 4.5. Assume, in addition, that equation f{z) = possesses at least two solutions vi 
and V2 such that f'{vi) > 0, z = 0, 1. Then, Vi will be exponentially stable equilibria of equation 
(|4.ip . Let now vq{x) be a smooth function such that 

Vo{x) = Vl, X £ Ql, Vo{x) = V2, X G 0,2 

for some non-empty fi, C 17 of the nonzero measure. 

Finally let us consider the initial data for problem (12. ip of the form 

Co := {vo{x),0,0). 

Then, since a > is small and the equilibria Vi are exponentially stable, the solution v{t,x) 
will remain close to Vi (for x € ilj) for all t. This shows that the smooth trajectory S{t)S,o = 
{v{t),dtv(t),w{t)) tends as t — > cxd to the discontinuous equilibrium (in the L^'-topology, accord- 
ing to the last theorem). 

This example shows that we cannot extend the assertion of the theorem to the case p = oo 
and obtain the convergence in the topology of the phase space <I>oo- Indeed, if the sequence 
of continuous functions converges in L°° to some limit function, this function is automatically 
continuous. Thus, the (j-limit set of the above constructed trajectory in the topology of the 
phase space is empty: 

Remark 4.6. It is clear from the proof of Proposition 14.21 that all of the equilibria TZa are 
hyperbolic in the phase space $oo- Thus, we may construct the infinite-dimensional stable and 
unstable manifolds for any equilibrium belonging to TZa if a > is small enough. However, 
it does not help much for the study of the limit dynamics since, as shown in the previous 
proposition, generically, we do not have the stabilization in the topology of ^oo, but only in a 
weaker space [LP(0)]^. And in this weaker space the solution semigroup S{t) is not differentiable. 
By this reason, we are not able to extract the exponential convergence from the hyperbolicity of 
any equilibrium and do not know whether or not such exponential convergence takes place. 

To conclude, we note that, arguing analogously to the proof of Proposition l4.2t one can extract 
some reasonable information about the equilibria TZa even in the case where a is not small. 

Proposition 4.7. Let the assumptions of Theorem \2. 6\ hold. Assume, in addition, that {v,iv) G 
is a spatially homogeneous hyperbolic (in <I>ooJ equilibrium of problem ()2.ip . i.e., that the 
equation 

(4.24) f'(^y)0-a{-A, + l)-^e = h 

is uniquely solvable for every h G L°°(il). Assume, finally, that there exists another constant 
V V such that 

(4.25) f{v) = f{v) 

(this, of course, may happen only in the case of non-monotone f). Then, there exists 5o > 
such that, for any measurable partition = $7i U $^2 on two disjoint sets where 

(4.26) |f^2|<<5o, 

there exists a hyperbolic equilibrium {v,w) such that v is close (in the L°° -metric) to 

(4.27) vi2 := vxQ.Ax) + vXn2{x)- 
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Proof. We first check that the equation of variations 

(4.28) f'{vi2)e-a{-A, + l)-^9 = h 

is uniquely solvable if the measure of is small. To this end, we construct the approximative 
solution of this equation in the form 9 := 9q + 9 where solves equation ()4.24p and 

(4.29) 9{x) := [f\v) - f'{d)]9o{x)xn,ix). 
Then, since u is a hyperbolic equilibrium, we have 

(4-30) ll^o||L°°(n) < C||^llL°°(n)- 

Moreover, the approximate solution 9 thus constructed solves 

(4.31) /'(^^i2)^ - a(-A, + iy^9 = h:= -a(-A, + l)-^9. 
Finally, since the measure of ^2 is small, we have 

||^||l2 < c|r22 1 "'"'^^ 11^0 — CiSq^ ||/i||/^tx) 

and, consequently, 

(4.32) ||/i||loo < C\\h\\H2 < C'\\9\\l2 < CaJ^^H/iIlL- , 

where the constant C2 is independent of h and of the concrete form of the partition Q = U$72- 

1/2 II ~ II 

Thus, if 5o > is so small that C26Q := k < 1, the norm of the reminder ||/i||loo is estimated 

k||/i||loo with K < 1. Then, the standard iteration process gives the desired solution 9 of equation 

(|4.28p together with estimate 

(4.33) II^IIl- < C||/i||lco 

with the constant C independent of 60 0. The uniqueness of a solution can be obtained in a 
standard way using the observation that the operator /'(f 12) — C({—Ax + 1)~^ is self-adjoint in 

It is now not difficult to finish the proof of the proposition. Indeed, we seek for the desired 
equilibrium {v,w) in the form 

v{x) = vuix) + 9{x), 
where is a small corrector which should satisfy the equation 

(4.34) f{vi2 + 9)- f{vi2) - a(-A, + l)-^9 = h := a(-A, + l)[{v - v)xn2] 
and, arguing as before, we see that 

< c\n2\^/^. 

Applying now the implicit function theorem to equation ()4.34p . we establish the existence of a 
unique solution 9, 

< Ci\n2\ ' 

measure of 0,2 is small enough. Proposition 14.71 is proved. □ 

Remark 4.8. Although we are not able neither to give a complete description of the equilibria 
set TZa nor to verify the stabilization if a is not small, we see that, under the assumptions of the 
last proposition (which are, in a sense, natural for the non-monotone case), the set of equilibria 
is not compact in $00 and is not closed in the weak-star topology of the phase space. These 
facts do not allow to extend the global attractor theory for the non- monotone case. 

Note also that, although we formulate (for simplicity) Proposition 14. 71 for the case of spatially- 
homogeneous hyperbolic equilibrium (v, w) it can be easily extended to the case of non-homo- 
geneous equilibria. This shows that the conclusion of Proposition W7l\ is somehow "generic" for 
the non-monotone case. 
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5. Appendix. Estimation of time derivatives for non-autonomous perturbations 

of regular attractors. 

The aim of that appendix is to verify the auxihary estimate for non-autonomous perturbations 
of regular attractors. To be more precise, consider an ODE in M": 

(5.1) u'{t) = F{u{t)), u{0) = no 

for some F G C^(M",M"). We assume that, for every uq G M", this equation is globahy (for 
t > 0) solvable and the associated semigroup S{t)uQ := u{t) is dissipative, i.e., 

(5.2) \\S{t)uo\\<Q{\\uo\\)e-^' + C, 

for some positive j3 and C=k and monotone Q. Therefore, equation (j5.ip possesses a global 
attractor A in M". Our main assumption is that this attractor is regular in the sense of Theorem 
13.81 i.e., all of the equilibria uq G TZ are hyperbolic, every trajectory, belonging to the attractor 
^ is a heteroclinic orbit connecting two different equilibria and the attractor ^ is a finite union 
of the unstable manifolds J^uo associated with the equilibria uq G TZ: 

(5.3) A = U.,,enMt,. 

Finally, we assume that the so-called no-cycle condition is satisfied, i.e., the attractor A does not 
contain any heteroclinic cycles. As known, that is always true in the case when (15.11) possesses 
a global Lyapunov function. 

Consider now the following small non-autonomous perturbation of equation ()5.ip 

(5.4) u' = F{u) + h{t), u(0) = no, t> 0, 
where the non-autonomous external force is uniformly small: 

(5.5) ll^lll4^i.°°(M+,R") < £ < 1- 
The main result of this appendix is the following estimate. 

Proposition 5.1. Let the above assumptions hold and let the external force h £ W^'°^{W^) 
satisfy estimate (jS.Sp for sufficiently small e > 0. Then, any solution u{t) of the perturbed 
problem (j5.4p satisfies the following estimate: 

(5.6) / \\u'{t)\\dt <Ci + C2 [ \\h'it)\\dt, 

Jo Jo 

where the positive constants Ci and C2 depend only on the norm of n(0) and are independent 

of T and the concrete choice of n(0) and h{t) . 

Proof. Indeed, using the standard regular attractor perturbation arguments, one can check that 
for every bounded set B of M" and every 5 > 0, there exist T = T{B,6) and eo = eo{B,6) 
such that, for every £ < Eq and every trajectory n(t) starting from B, we can find a sequence 
n^ • • • , n^v of different equilibria Ui € TZ (of problem (15. ip !) and a sequence of times 

= T(+ < Tf < T+ < < < • • • < < r+ = 00 

such that 

(5.7) u{t)eOs{ui), te{Tr,T+), T- - < T, i = I, ■ ■ ■ , N. 

In other words, the sequence of equilibria Ui and the values of Tj^ depend on the concrete choice 
of n(0) and h, but the number N of equilibria is bounded by the whole number ^TZ < 00 of 
possible equilibria (since the equilibria must be different) and the differences — T^" ^ are also 
uniformly bounded by T, see O El [20] for the details. 

Thus, any trajectory starting from B spends at most time Tout ■= H^TZ ■ T outside of the 
5- neighbourhood ©^(M) of the equilibria set TZ and this time depends only on B and 6. By 
this reason, the part of the trajectory, lying outside of ©^(M) gives only a finite and uniformly 
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bounded impact to the integral (j5.6p (which can be included into the constant Ci). So, we 
only need to estimate the left-hand side of (j5.6p for the case where u{t) belongs to a small 
neighbourhood of a single fixed equilibrium uq & TZ only. 

To this end, we will use the hyperbolicity assumption on uq. Indeed, the implicit function 
theorem implies the existence of > and 6 > such that, for every e < eo> there exists a 
unique solution UuQ,h{t) of (15. 4p belonging to the (5-neighbourhood of uq for all t. Moreover, this 
solution, in a fact, belongs to the Ce- neighbourhood of uq and the following estimate holds: 

(5.8) \K,,dt)\<C [ e'^\'-^\\h'is)\ds, 

where the constant C and the hyperbolicity exponent k are independent of the concrete choice 
of Uq £ TZ and the external force h satisfying (|5.5p , see [TUl [20] . 

Furthermore, since uq is hyperbolic, the trajectory UuQ,h{t) will be also hyperbolic and we 
will have an exponential dichotomy in a small (^-neighbourhood of UuQ,h- In particular, every 
trajectory u{t) belonging to Os{uq) for t G [0,5], S ^ 1, will tend exponentially to UuQ,h{t) 
inside of the interval 

(5.9) \u{t) - Uu„H{t)\ + \n{t) - K^^dt)\ < C(e-«l*l + e-'^l^-*!) 

and C and k are independent of the concrete choice of u and h, see [TOl ED] for the details. 
Therefore, 

(5.10) f\u'{s)\ds<C+ f\K,,^{s)\ds 

Jo Jo 

for t G [0, S] and u{t) £ Os{uo). 
Thus, we have proved that 

(5.11) r \u'is)\ ds<ci+c2Y. r \K„h{^)\ ds, 

where the constants Cj depend only on the radius of B. 

In order to deduce (15.61) from (j5.1ip . we will use estimate (15. Sp . Indeed, integrating it over 
t £ [0,T] and using that h can be extended for t < by zero, we have 

(5.12) / \K^h{t)\dt<C \h'{t)\dt + C e-''\^-^\\h'{s)\ds <Ci + C \h'{t)\dt. 
Jo ' Jo Jt Jo 

Inserting this estimate into the right-hand side of (j5.1ip , we obtain (|5.6p and finish the proof of 
the proposition. □ 

Remark 5.2. As we can see from the proof, estimate (j5.6p has a general nature which is not 
restricted by the class of ordinary differential equations. However, since we use it in the paper 
only for the ODE case, in order to avoid the technicalities related with the formulation of a 
"general" PDE, we restrict our consideration to the case of a ODE as well. 
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